numerical simulation of plume rise in neutral atmospheric stability condition

نویسندگان

خسرو اشرفی

استادیار گروه مهندسی عمران محیط زیست دانشکدۀ محیط زیست دانشگاه تهران مجید شفیع پور مطلق

استادیار گروه مهندسی عمران محیط زیست دانشکدۀ محیط زیست دانشگاه تهران حدیث عباس زاده امیردهی

کارشناس ارشد مهندسی عمران محیط زیست دانشکدۀ محیط زیست دانشگاه تهران

چکیده

introduction air pollution has harmful effect on human health and the environment. accordingly, considerable effort has been put to analyze the air pollutants. one important issue is the spatial distribution of these pollutants. dispersion of the pollutants released from sources on the ground is mostly driven by the planetary boundary layer where turbulent flow causes mixing of the content of the stationary ground layer of atmosphere with higher moving layers and thereby clears out the pollutants rapidly. to predict the dispersion of air pollution in the atmosphere, researchers investigate the behavior of plume rise under certain conditions. two factors critically affect the result of the test: primary plume rise and its dispersion. in this work, we study the effective height of emission which determines the dispersion of the pollutants. the height of the plume has been subjected to semi-empirical and to more accurate numerical studies. these approaches suffer from some shortcomings. for instance, the gravitational effects and fluctuations of the wind speed are neglected. numerical approaches are typically based on a gaussian plume rise model so called isc3 (published by epa). in this paper plume rise of emission of an air pollutant is simulated using fluent software which allows one to input natural wind velocity profile and to consider the gravity. another advantage of this approach lies in the usage of turbulent kinetic energy (tke) and temperature profiles. in this work, theoretical plume rise using the isc3 model were being calculated at first. the calculated values of xf and ∆h are equal to 437.85 m and 50.87m, respectively. then, the dispersion was being simulated accordingly using fluent and the results were compared with those of the numerical studies.      materials and methods to simulate the plume rise, first a sample stack has been designed and meshed in gambit and then transferred to fluent.without loss of the generality of the approach, the following assumptions were made for physical factors of the chimney and weather condition. the simulation was performed in two dimensions (2d) on an area with a length of 2000 m and a height of 1000 m. the opening of the chimney was modeled at a height of 60 m and in an area of 500 m far from the entrance. the exit velocity is assumed to be 10m/s with a temperature of 450 k. the wind velocity at 10 m height is 4 m/s. generating a good mesh is an essential part of any cfd problem. we generated the geometry by the use of gambit, the preprocessing module of the fluent. in the present simulation, a standard k-ε model is used and neutral stability condition of the atmosphere is setup using user defined function (udf) facility of the fluent software by defining temperature, wind, and turbulence kinetic energy (tke) profiles at the inlet of the atmospheric boundary layer. the wind velocity profile was obtained from aermod and has a logarithmic form. the vertical profile of the temperature for neutral stability condition is as follow:         t= -0.0098 δz + t0                            (1) where, and represent the temperature and the elevation gradient, respectively. one of the significant steps is to insert the turbulence kinetic energy (tke) profile which depends on heat transfer, momentum, and moisture in the boundary layer and the tke obtained from the following try and error method. in order to evaluate the effect of the gravitational acceleration on the plume rise, the simulation has been done with considering the gravitational acceleration and then the results have been compared to the results when the gravitational acceleration has been ignored.   results and discussions to better compare the simulation results with the ones for semi-empirical equations, the maximum concentration and the geometric mean concentration have been averaged. the plume rise with the effect of the gravitational acceleration is shown in figure 1. simulation repeated with ignoring the effect of the gravitational acceleration and the results were shown in figure 2. the maximum concentration reached its highest value approximately at 20 meter far from the stack with an effective plume height of 2.43 m when the gravitational acceleration affected the flow of the plume. but as shown in figure 2, when gravitational acceleration has been ignored in the simulation, the plume height reached to its highest point at about 450 meter far from the stack with an effective plume height of 8.43 m.     figure 1. plume rise in neutral atmospheric stability with the effect of gravitational acceleration   figure 2. plume rise in neutral atmospheric stability with ignoring the effect of gravitational acceleration     conclusions the quantity of ∆has a function of xf  was presented in table 1. ∆h is obtained from the average of maximum concentration and geometric mean concentration. for the neutral atmospheric stability condition, the results give better agreement with the semi-empirical equations when the gravitational acceleration has been ignored. but, the plume rise became 2.5 times larger than the values obtained from the semi-empirical equations when the gravitational acceleration was applied.  this occurs because the gravitational acceleration drawn most mass of the outlet flow; therefore, buoyancy force applied to less concentration of the outlet flow and dispersed it much more in the vertical and horizontal directions. as the presented simulation considers actual wind velocity profile, actual temperature profile, tke profile, and gravitational acceleration effect, the results are more close to the reality. table 1. comparison of different approaches for estimation of plume rise plume rise (with the effect of gravitational acceleration) plume rise (ignoring the effect of gravitational acceleration) semi-empirical relations   128.94 84.53 50.87 h∆ 437.85 437.85 437.85 xf

برای دانلود باید عضویت طلایی داشته باشید

برای دانلود متن کامل این مقاله و بیش از 32 میلیون مقاله دیگر ابتدا ثبت نام کنید

اگر عضو سایت هستید لطفا وارد حساب کاربری خود شوید

منابع مشابه

Boundary Condition in Particle-based Visual Simulation of Plume Fluid

Recently, in the field of CG, the animation of flame or smoke which uses the fluid simulation is researched. In this paper, we researched the particle-based visual simulation of plume fluid. Plume is fluid with a mushroom-like shape. For example, it is explosive flame or cumulus. Plume, such as continuously growing flame or cloud with an expansive gush, is a difficult phenomenon to represent wi...

متن کامل

Sensitivity of air quality simulation to smoke plume rise

Plume rise is the height smoke plumes can reach. This information is needed by air quality models such as the Community Multiscale Air Quality (CMAQ) model to simulate physical and chemical processes of point-source fire emissions. This study seeks to understand the importance of plume rise to CMAQ air quality simulation of prescribed burning to plume rise. CMAQ simulations are compared between...

متن کامل

Wellbore Stability in Shale Formation Using Analytical and Numerical Simulation

Optimization of drilling fluid parameters such as mud weight, salt concentration, and temperature is essential to alleviate instability problems during drilling through shale sections. The selection of suitable mud parameters can benefit from analyses that consider significant instability processes involved in shale-drilling–fluid interactions. This paper describes the development of analytical...

متن کامل

Including the sub-grid scale plume rise of vegetation fires in low resolution atmospheric transport models

1 Center for Weather Forecasting and Climate Studies, INPE, Cachoeira Paulista, Brazil 2 NASA Ames Research Center, Moffet Field, USA 3 USDA Forest Service, Montana, USA 4 Max Planck Institute for Chemistry, Mainz, Germany 5 UW-Madison Cooperative Institute for Meteorological Satellite Studies, Madison, WI, USA 6 Laboratório de Combustão e Propulsão, INPE, Cachoeira Paulista, Brazil 7 Departmen...

متن کامل

Numerical Simulation of the Cathode Plume of a Hall Thruster

The plasma plume generated by the cathode of a 6-kW laboratory Hall thruster is modeled using a hybrid particle-fluid method. The approach utilizes the direct simulation Monte Carlo method to simulate the collision dynamics of heavy particles, the particle-in-cell method to model the transport of heavy species including the effects of electrostatic fields on the ions, and a “fluid” method that ...

متن کامل

Wildland fire emissions, carbon, and climate: Plume rise, atmospheric transport, and chemistry processes

This paper provides an overview and summary of the current state of knowledge regarding critical atmospheric processes that affect the distribution and concentrations of greenhouse gases and aerosols emitted from wildland fires or produced through subsequent chemical reactions in the atmosphere. These critical atmospheric processes include the dynamics of plume rise, chemical reactions involvin...

متن کامل

منابع من

با ذخیره ی این منبع در منابع من، دسترسی به آن را برای استفاده های بعدی آسان تر کنید


عنوان ژورنال:
محیط شناسی

جلد ۴۰، شماره ۱، صفحات ۲۳۷-۲۴۶

کلمات کلیدی
introduction air pollution has harmful effect on human health and the environment. accordingly considerable effort has been put to analyze the air pollutants. one important issue is the spatial distribution of these pollutants. dispersion of the pollutants released from sources on the ground is mostly driven by the planetary boundary layer where turbulent flow causes mixing of the content of the stationary ground layer of atmosphere with higher moving layers and thereby clears out the pollutants rapidly. to predict the dispersion of air pollution in the atmosphere researchers investigate the behavior of plume rise under certain conditions. two factors critically affect the result of the test: primary plume rise and its dispersion. in this work we study the effective height of emission which determines the dispersion of the pollutants. the height of the plume has been subjected to semi empirical and to more accurate numerical studies. these approaches suffer from some shortcomings. for instance the gravitational effects and fluctuations of the wind speed are neglected. numerical approaches are typically based on a gaussian plume rise model so called isc3 (published by epa). in this paper plume rise of emission of an air pollutant is simulated using fluent software which allows one to input natural wind velocity profile and to consider the gravity. another advantage of this approach lies in the usage of turbulent kinetic energy (tke) and temperature profiles. in this work theoretical plume rise using the isc3 model were being calculated at first. the calculated values of xf and ∆h are equal to 437.85 m and 50.87m respectively. then the dispersion was being simulated accordingly using fluent and the results were compared with those of the numerical studies.      materials and methods to simulate the plume rise first a sample stack has been designed and meshed in gambit and then transferred to fluent.without loss of the generality of the approach the following assumptions were made for physical factors of the chimney and weather condition. the simulation was performed in two dimensions (2d) on an area with a length of 2000 m and a height of 1000 m. the opening of the chimney was modeled at a height of 60 m and in an area of 500 m far from the entrance. the exit velocity is assumed to be 10m/s with a temperature of 450 k. the wind velocity at 10 m height is 4 m/s. generating a good mesh is an essential part of any cfd problem. we generated the geometry by the use of gambit the preprocessing module of the fluent. in the present simulation a standard k ε model is used and neutral stability condition of the atmosphere is setup using user defined function (udf) facility of the fluent software by defining temperature wind and turbulence kinetic energy (tke) profiles at the inlet of the atmospheric boundary layer. the wind velocity profile was obtained from aermod and has a logarithmic form. the vertical profile of the temperature for neutral stability condition is as follow:         t= 0.0098 δz + t0                            (1) where and represent the temperature and the elevation gradient respectively. one of the significant steps is to insert the turbulence kinetic energy (tke) profile which depends on heat transfer momentum and moisture in the boundary layer and the tke obtained from the following try and error method. in order to evaluate the effect of the gravitational acceleration on the plume rise the simulation has been done with considering the gravitational acceleration and then the results have been compared to the results when the gravitational acceleration has been ignored.   results and discussions to better compare the simulation results with the ones for semi empirical equations the maximum concentration and the geometric mean concentration have been averaged. the plume rise with the effect of the gravitational acceleration is shown in figure 1. simulation repeated with ignoring the effect of the gravitational acceleration and the results were shown in figure 2. the maximum concentration reached its highest value approximately at 20 meter far from the stack with an effective plume height of 2.43 m when the gravitational acceleration affected the flow of the plume. but as shown in figure 2 when gravitational acceleration has been ignored in the simulation the plume height reached to its highest point at about 450 meter far from the stack with an effective plume height of 8.43 m.     figure 1. plume rise in neutral atmospheric stability with the effect of gravitational acceleration   figure 2. plume rise in neutral atmospheric stability with ignoring the effect of gravitational acceleration     conclusions the quantity of ∆has a function of xf  was presented in table 1. ∆h is obtained from the average of maximum concentration and geometric mean concentration. for the neutral atmospheric stability condition the results give better agreement with the semi empirical equations when the gravitational acceleration has been ignored. but the plume rise became 2.5 times larger than the values obtained from the semi empirical equations when the gravitational acceleration was applied.  this occurs because the gravitational acceleration drawn most mass of the outlet flow; therefore buoyancy force applied to less concentration of the outlet flow and dispersed it much more in the vertical and horizontal directions. as the presented simulation considers actual wind velocity profile actual temperature profile tke profile and gravitational acceleration effect the results are more close to the reality. table 1. comparison of different approaches for estimation of plume rise plume rise (with the effect of gravitational acceleration) plume rise (ignoring the effect of gravitational acceleration) semi empirical relations   128.94 84.53 50.87 h∆ 437.85 437.85 437.85 xf

میزبانی شده توسط پلتفرم ابری doprax.com

copyright © 2015-2023